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This  study  investigates  a  robust  controller  in  regulating  the  pulse  width  modulation  (PWM)  of  a  DC/DC 
converter  for  a  polymer  electrolyte  membrane  fuel  cell  (PEMFC)  application.  A  significant  variation  in  the 
output  voltage  of  a  PEMFC  depends  on  the  power  requirement  and  prevents  a  PEMFC  from  directly 
connecting  to  a  subsequent  power  bus.  DC/DC  converters  are  utilized  to  step-up  or  step-down  voltage  to 
match  the  subsequent  power  bus  voltage.  In  this  study,  a  full  dynamic  model,  which  includes  a  PEMFC 
and  boost  and  buck  DC/DC  converters,  is  developed  under  MATLAB/Simulink  environment  for  control.  A 
robust  PWM  duty  ratio  control  for  the  converters  is  designed  using  time  delay  control  (TDC).  This  control 
enables  state  variables  to  accurately  follow  the  dynamics  of  a  reference  model  using  time-delayed  in¬ 
formation  of  plant  input  and  output  information  within  a  few  sampling  periods.  To  prove  the  superiority 
of  the  TDC  performance,  traditional  proportional-integral  control  (PIC)  and  model  predictive  control 
(MPC)  are  designed  and  implemented,  and  the  simulation  results  are  compared.  The  efficacies  of  TDC  for 
the  PEMFC-fed  PWM  DC/DC  converters  are  validated  through  experimental  test  results  using  a  100  W 
PEMFC  as  well  as  boost  and  buck  DC/DC  converters. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Fuel  cell  is  a  potential  clean  energy  solution  to  remedy  the 
increasing  environmental  pollution  issues,  and  the  polymer  elec¬ 
trolyte  membrane  fuel  cell  (PEMFC)  shows  superiority  in  the 
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automobile  industry  [1—3].  The  PEMFC  exhibits  higher  efficiency, 
lower  or  zero  emission,  lower  noise,  and  lower  operation  temper¬ 
ature  compared  with  other  types  of  fuel  cells.  For  vehicular  appli¬ 
cation,  a  regulated  output  voltage  of  PEMFC  is  compulsory  to  feed 
the  power  for  the  subsequent  electrical  motor  system;  however,  the 
output  voltage  of  a  PEMFC  is  not  constant  because  of  the  charac¬ 
teristic  of  its  polarization  curve,  which  is  governed  by  an  electro¬ 
chemical  reaction.  The  electrochemical  reaction  is  affected  by 
variations  in  the  operating  points  of  parameters  such  as  mass  flow 
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rate,  inlet  pressure,  stack  temperature,  water  content,  relative  hu¬ 
midity,  stoichiometry  of  air  and  hydrogen,  and  channel  geometry 
[4-10].  Moreover,  recent  power  buses  for  fuel  cell  hybrid  vehicles 
require  different  dual  voltage  links,  namely,  higher  and  lower  links. 
A  higher  link  of  24  V  is  used  for  electric  motors,  air-conditioners, 
and  pump  connections,  whereas  the  lower  link  of  12  V  is  used  for 
conventional  switches  and  lamp  connections  [3].  Therefore,  DC/DC 
converters  are  the  best  choice  for  regulating  the  output  voltage  of 
PEMFC  to  satisfy  this  dual  power  bus  requirement.  Additionally,  DC/ 
DC  converters  can  step-up  or  step-down  input  voltage  depending 
on  the  different  operating  conditions,  i.e.,  boost  converter  for  high 
voltage  (step-up  voltage)  [8,11,12]  and  buck  converter  for  low 
voltage  (step-down  voltage)  [13,14].  Thus,  PEMFC-fed  DC/DC  con¬ 
verter  system  can  be  regarded  as  a  typical  integrated  power  supply 
system  (Fig.  1)  that  can  constantly  supply  the  desired  DC  link 
power. 

Controlling  the  integrated  power  system  exhibits  challenging 
problems  because  PEMFC  and  DC/DC  converters  are  non-linear 
systems  resulting  in  coupled  characteristics.  To  control  the  inte¬ 
grated  system,  the  construction  of  a  simple  and  accurate  control- 
oriented  dynamic  model  is  necessary.  Choe  et  al.  [8]  developed 
an  associated  PEMFC  pulse  width  modulation  (PWM)  DC/DC  con¬ 
verter  model  by  considering  balance-of-plant  (BOP)  of  a  PEMFC, 
including  air  supply  and  thermal  systems,  and  a  DC/DC  boost 
converter.  This  model  showed  comprehensive  control  variables 
with  high-order  dynamics  and  a  complex  control  process  validated 
only  by  simulation  results.  A  more  detailed  PEMFC  model  for 
employing  DC/DC  converter  can  be  found  in  Refs.  [15,39].  To 
simplify  the  modeling  process  and  reduce  computational  burden 
for  real-time  control  implementation,  Bjazic  et  al.  11  ]  and  Shtessel 
et  al.  [22]  presented  equivalent  electrical  circuit  dynamic  models  of 
a  PEMFC,  which  were  subsequently  used  in  a  DC/DC  boost  con¬ 
verter  with  a  double  layer  capacitor  to  simulate  the  electrode 
behavior  of  a  PEMFC  system  16].  A  modified  equivalent  electrical 
circuit  approach  with  variable  capacitance  was  implemented  in 
Ref.  [12]. 

An  accurate  and  robust  PWM  duty  ratio  control  is  required  to 
manipulate  output  voltage  that  precisely  tracks  the  setpoint 
value,  because  DC/DC  converters  are  inherently  non-linear  with 
variations  in  system  parameters.  Numerous  studies  have  reported 
on  the  accurate  and  robust  PWM  duty  ratio  control  for  DC/DC 
converter.  Traditional  small  signal  state-space  averaging  method 
was  first  summarized  in  Ref.  [17].  This  method  is  widely  used  for 
the  linearization  power  converter  to  construct  linear  feedback 
control  systems,  such  as  linear  quadratic  regulator  control  [18] 
and  proportional-integral  (PI)  or  proportional-integral- 
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Fig.  1.  PEMFC  integrated  power  supply  system. 


derivative  control  [8,12,19].  However,  the  linear  control  of  DC/ 
DC  converter  is  not  suitable  for  PEMFC  application  because  the 
PEMFC  operation  point  is  always  varying  compared  with  the 
normal  DC  power  source.  Many  robust  control  approaches,  such 
as  Hoc  control  [20],  gain-scheduled  control  [21],  and  sliding  mode 
control  (SMC),  have  recently  been  developed  to  cope  with  non¬ 
linearity  and  ensure  stability  [13,22-25].  Moreover,  several 
complicated  control  techniques  have  been  introduced  to  improve 
the  DC/DC  converter  performance.  Wai  et  al.  [26]  presented  a 
total  SMC  with  adaptive  linear  neural  network  filter  to  reduce 
current  ripple  to  protect  clean-energy  sources.  Geyer  et  al.  [27] 
introduced  a  model  predictive  control  (MPC)  based  on  hybrid 
piecewise  affine  model  approximations  to  efficiently  reject  dis¬ 
turbances  inherent  in  the  input  voltage.  Similar  results  were 
collected  in  Ref.  [28],  in  which  five  techniques  from  hybrid  and 
optimal  control  were  evaluated  and  compared  to  diminish  the 
converter  output  voltage  ripples  under  variations  in  supply 
voltage  and  load  resistor.  The  dynamic  evolution  control  for  fuel 
cell-based  DC/DC  converter  that  satisfies  the  desired  performance 
was  developed  in  Ref.  [29].  Optimization  algorithms  such  as 
neural  network  (NN)  [30  and  genetic  algorithm  31]  have  also 
been  reported  to  tune  the  control  parameters  to  improve  DC/DC 
converter  control  performance. 

This  study  presents  a  new  design  methodology  of  time  delay 
control  (TDC)  for  the  PEMFC-fed  PWM  boost  and  buck  DC/DC 
converter  control.  TDC  is  a  totally  different  control  algorithm  from 
the  widely  used  robust  control  concepts,  such  as  SMC  and  NN, 
because  it  is  based  on  the  model  reference  control.  TDC  was  orig¬ 
inally  proposed  by  Youcef-Toumi  et  al.  [32].  The  prominent 
advantage  of  TDC  is  its  ability  to  represent  unknown  system  dy¬ 
namics  or  disturbance  using  the  time-delayed  information  of  plant 
input  and  output  relation  just  within  several  sampling  periods.  TDC 
control  system  structure,  stability,  controllability,  and  observability 
have  been  discussed  in  some  studies  [33,34].  Moreover,  TDC  has 
been  widely  applied  in  various  industrial  areas,  such  as  robot 
manipulation  control  [35],  PEMFC  air  supply  system  implementa¬ 
tion  [15,36],  and  boost  converter  control  [37]. 

In  the  present  study,  an  integrated  system  dynamic  model  is 
derived  in  Section  2,  where  an  equivalent  electrical  circuit  with 
variable  double-layer  capacitor  is  introduced  to  replace  PEMFC 
dynamic  characteristic.  TDC  is  utilized  to  design  the  DC/DC  con¬ 
verter  PWM  duty  ratio  controller  in  Section  3.  In  Section  3,  two 
conventional  control  methods,  namely,  PIC  and  MPC,  are  performed 
on  the  same  platform  to  show  the  superiority  of  TDC.  Simulation 
results  of  the  designed  PEMFC-fed  PWM  DC/DC  converters  are 
expressed  in  Section  4.  To  validate  the  efficacy  of  TDC,  experiments 
based  on  boost  and  buck  DC/DC  converters  with  commercial 
PEMFC  Horizon  H-100  with  National  Instrument  (NI)  DAQPCI  board 
for  real-time  control  are  performed,  and  the  results  are  also  pre¬ 
sented  in  this  section. 

2.  PEMFC-based  DC/DC  converter  model 

PEMFC-based  DC/DC  converter  model  requires  the  internal 
relation  between  PEMFC  and  DC/DC  converter  because  the  output 
voltage  of  PEMFC  is  the  supply  voltage  of  the  DC/DC  converter. 
Therefore,  the  development  of  a  full  control-oriented  model 
including  PEMFC  and  DC/DC  converters  is  necessary.  The  highest 
accuracy  of  PEMFC  model  can  be  obtained  from  a  CFD-based  model 
with  numerous  state  variables  [6,10]  or  from  a  comprehensive  BOP- 
based  model  [15].  However,  these  approaches  cannot  be  fully 
managed  in  practice  because  of  the  large  amount  of  highly  coupled 
variables.  Therefore,  these  approaches  are  not  suitable  for  control. 
The  major  characteristic  of  PEMFC  is  indicated  by  its  polarization 
curve,  which  intuitively  shows  the  PEMFC  output  voltage  versus 
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stack  current  density  relation.  Therefore,  equivalent  electrical  cir¬ 
cuit  models  with  variable  double-layer  capacitor  can  be  used  in  our 
study  to  control  the  PEMFC  electrochemical  transient  characteris¬ 
tics.  This  model  can  facilitate  controller  design  with  high  accuracy 
and  small  number  of  state  variable.  The  validation  of  a  developed 
PEMFC  model  is  proven  by  comparing  the  test  results  obtained 
from  Horizon  H-100,  a  commercial  PEMFC,  which  comprised 
20  cells  with  100  W  rated  power  [38].  H-100  is  a  convenient  and 
delicate  PEMFC  that  is  air  blown  and  cooled  by  an  integrated  fan, 
supplied  with  hydrogen,  purged  by  two  valves,  and  equipped  with 
self-humidification.  Thus,  it  is  suitable  for  low-power  system 
application. 


2  A.  PEMFC  model 
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where  l  and  A  are  the  membrane  thickness  and  cell  active  area, 
respectively.  Am  is  the  average  water  content  in  the  membrane, 
which  is  based  from  the  membrane  humidity  [8]. 

The  last  term  of  Eq.  (1)  corresponds  to  the  concentration  loss 
Vcon,  caused  by  the  insufficient  reactant  delivery  when  the  current 
density  occurs  and  is  given  by  [16] 

v”"  -  ->{'  -£)  (8> 

where  j  is  fuel  cell  current  density  defined  by 


The  principle  of  PEMFC  reaction  involves  the  transport  of 
hydrogen  to  the  anode  catalyst  layer  and  its  split  into  protons  and 
electrons  under  a  defined  pressure.  The  membrane  can  transfer 
protons  only  to  the  cathode,  but  it  can  stop  electrons.  Subsequently, 
electrons  generate  electricity  through  an  external  electrical  circuit. 
Protons,  oxygen,  and  electrons  react  to  produce  water  and  heat  at 
the  cathode.  The  output  voltage  of  the  single  cell  is  described  by  [9] 

^cell  =  ^Nernst  —  ^act  —  ^ohm  —  ^con  ( 1 ) 

where  ENernst  is  the  reversible  voltage,  Vact  is  the  activation  drop, 
Vohm  is  the  ohmic  drop,  and  VCon  is  the  concentration  drop.  The 
reversible  voltage  is  calculated  by  Nernst  potential  equation  [13].  Its 
value  depends  on  operation  temperature  and  pressure.  For  a 
standard  condition  (T0  =  298.15  K,  P0  =  101,325  Pa),  the  standard 
reversible  voltage  is  1.229  V  [15].  Then,  the  Nernst  potential 
equation  can  be  written  as 


The  electrochemical  transient  characteristics  of  PEMFC  can  be 
simulated  by  an  equivalent  electrical  circuit,  which  is  based  on  the 
double-layer  capacitor  charging  effect.  This  effect  results  in  the 
output  voltage  with  a  small  time  lag  when  the  current  suddenly 
changes.  This  phenomenon  is  precisely  represented  using  the  first- 
order  time  constant  made  by  the  capacitor  in  the  activation  and 
concentration  polarization  [8,16].  The  modified  equivalent  circuit 
with  variable  capacitor  model  [13]  is  utilized  in  this  study  because 
it  exhibits  better  experimental  fitting  performance.  Fig.  2  shows  the 
equivalent  electrical  circuit  of  a  PEMFC  dynamic  model,  and  the 
equivalent  capacitance  is  determined  by 
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Based  on  the  simplified  parametric  equivalent  of  the  general 
Butler- Volmer  equation,  the  activation  drop  is  given  by  [12,13] 


where  f  5  is  a  constant  value.  The  activation  polarization  resistance 
Pact  is  estimated  by 


d  _  ^act 
Kact  —  -■ — 
l  FC 


(ii) 


Subsequently,  the  voltage  of  the  capacitor  Vc  is  given  by 


Vact  =  -[?i+?2-T  +  ?3-Tln(Co2)+?4-nn(iFc)]  0) 

where  i>c  is  the  fuel  cell  current,  and  flf  ?2.  ?3>  and  ?4  are  the  cell 
parametric  coefficients,  the  values  of  which  are  identified  by 
experimental  fitting.  Cq2  represents  the  concentration  of  oxygen 
dissolved  at  the  cathode  catalytic  interface  and  is  calculated  by 

P°2 

r  =  _ 101,325 _  (4) 

°2  5.08  X  106  exp (—498)  1  ; 

The  ohmic  loss  caused  by  the  resistance  of  proton  conduction 
through  the  membrane  and  electron  transfer  through  the  elec¬ 
trodes  is  given  by  [8,9] 

^ohm  =  *’fc(Km+Kc)  (5) 

where  the  equivalent  resistance  of  electron  transfer  Pc,  is  usually 
considered  as  a  constant  value.  Pm  represents  the  resistance  of 
proton  conduction  obtained  by 
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dVc 

dt 


(12) 


where  pu  is  the  specific  resistivity  denoted  by  [9] 


Fig.  2.  Equivalent  electrical  circuit  of  PEMFC  dynamic  model. 
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Accordingly,  the  output  voltage  of  the  single  cell  based  on  the 
simplified  dynamic  model  where  Vc,  instead  of  Vact,  is  obtained 
by 

pcell  —  ^Nernst  —  Vq  ~  ^ohm  —  ^con  (13) 

The  overall  stack  formed  by  N  cells  output  voltage  is  calculated 
by 

pfc  =  N  •  vce\\  (14) 

The  parameters  and  coefficients  determined  from  experimental 
fittings  are  listed  in  Table  1.  Fig.  3(a)  shows  the  polarization  and 
power  curve  comparisons  between  the  proposed  model  and 
experimental  results.  The  experimental  results  were  obtained 
through  200  tests  denoted  by  error  bars.  The  transient  responses  of 
the  stack  output  voltage  with  the  stack  current  input  variations  are 
also  indicated  in  Fig.  3(b).  The  fitting  between  the  model  and  the 
real  cell  shows  acceptable  accuracy.  Therefore,  it  can  be  used  as  a 
PEMFC  control-oriented  model  with  a  small  degree  of  freedom. 
Some  deviations  occurred  at  low-current  operation  condition 
during  the  transient  operation,  but  low  power  is  not  usually  used  in 
actual  PEMFC  operation. 


2.2.  DC/DC  converter  model 

The  schematics  of  DC/DC  converters  are  shown  in  Fig.  4,  which 
includes  boost-  and  buck-type  converters  with  PEMFC  equivalent 
circuit.  Each  converter  comprised  a  MOSFET  switch  S,  a  diode  D, 
an  inductor  L,  and  a  capacitor  C.  The  external  load  variable  is 
represented  by  r0.  After  converter  conditioning,  the  PEMFC 
output  voltage  vFc  can  be  regulated  as  the  desired  load-fed 
voltage  v0.  An  inductance  large  enough  to  preserve  sufficient 
energy  to  charge  or  discharge  the  output  capacitor  power  is 
assumed.  Hence,  the  discontinuous  conduction  operation  mode 
can  be  neglected. 

The  state-space  form  of  DC/DC  converter  is  described  as 

x  =  Ax  +  Bu  +  Gw  nc.x 

y  =  Cx  +  Du  1  j 

where  u  represents  the  MOSFET  switch  location  controlled  through 
PWM  duty  ratio,  taking  values  from  the  set  of  {0:1}.  y  is  the  output 
of  the  system  physically  calling  the  load  terminal  voltage  v0.  w  is  the 
supply  voltage  vFc  regarded  as  the  system  input.  A ,  B,  G,  C,  and  D  are 
the  system  dynamic  matrices  that  govern  converter  characteristics. 
For  boost  DC/DC  converter 


Table  1 

List  of  PEMFC  parameters  used  in  simulation. 


Parameter 

Parameter  description 

Value 

T 

Stack  temperature 

333.15  K 

ph2 

Partial  pressure  of  hydrogen 

2  x  101,325  Pa 

po2 

Partial  pressure  of  oxygen 

101325  Pa 

£i.  ?2,  h.  U  ;5 

Cell  parametric  coefficients 

-1.103,  3.48  x  10"3, 

5.8  x  10-5,  -9  x  10"5,  20 

Rc 

Constant  resistance 

2.05  x  10"2  Q 

l 

Thickness  of  the  membrane 

0.0178  cm 

A 

Area  of  the  membrane 

22.5  cm2 

Membrane  water  content 

11 

R 

Gas  constant 

8.3145  J  (mol  K)"1 

F 

Faraday  constant 

96,485  C  mol"1 

Jmax 

Maximum  current  density 

0.622  A  cm-2 

N 

Number  of  cell 

24 
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For  buck  DC/DC  converter 
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(17) 


Obviously,  controlling  a  boost  converter  is  more  difficult  than 
controlling  a  buck  converter.  Moreover,  the  boost  type  displays  the 
non-minimum  phase  phenomenon  with  the  output  voltage  and 
current  having  different  phase  characteristics  [28]. 


3.  PWM  control  approaches 

In  this  section,  three  different  control  approaches  are  proposed 
to  control  the  PWM  duty  ratio.  Tables  2  and  3  provide  the  circuit 
parameters  of  the  prototype  boost  and  buck  converters,  respec¬ 
tively.  TDC  is  first  mentioned,  and  two  conventional  approaches, 
namely,  PIC  and  MPC,  are  also  introduced  to  compare  the  control 
performances. 


3.1.  New  design  based  on  TDC 

A  general  form  of  a  non-linear  time  invariant  system  is 
described  by 

x(t)  =  f(x(t ),  t)  +  h(x(t),  t)  +  B(x(t),  t)u(t)  +  d(t)  (18) 

where  x(t)  is  the  state  vector,  f[x(t),t)  is  the  known  dynamics  vector, 
h(x{t),t)  is  the  unknown  dynamics  vector,  B(x(t),t)  is  the  input 
matrix,  u(t)  is  the  control  vector,  and  d(t)  is  the  disturbance  vector. 
Assuming  that  the  state  vectors  and  their  derivatives  (Eq.  (18)),  are 
all  observable,  the  tracking  performance  satisfies  the  reference 
model  defined  as 

Xm(t)  =  AmXm(t)  +  Bmr(t)  (19) 

where  xm{t)  and  r(t)  represent  the  reference  model  state  vector  and 
input  vector,  respectively.  Am  and  Bm  are  the  system  matrices. 
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Fig.  3.  Comparison  between  modeled  and  experimental  results:  (a)  PEMFC  polarization  and  power  curves;  (b)  Transient  responses  of  PEMFC  output  voltage. 


Defining  the  error  as  e(t)  =  xm(t )  -  x(t),  the  error  dynamics  is 
given  by  [32] 


Subsequently,  the  state  space  model  can  be  written  by 
substituting  ESR  =  0  into  Eqs.  (16)  and  (17): 


e(t)  =  Aee(t)  =  (Am  +  K)e(t)  (20) 

where  K  represents  the  error  feedback  gain  matrix.  The  error  sys¬ 
tem  is  asymptotically  stable  when  the  eigen  values  of  Ae  are  located 
at  the  left  s-  half  plane.  Combining  Eqs.  (19)  and  (20)  and  rear¬ 
ranging  the  terms  will  obtain 

B(x(t),  t)u(t)  =  -f(x(t),  t)  -  h(x(t),  t)  -  d(t )  +  xm(t)  -  Aee(t ) 

(21) 

If  B(x(t),t)  is  not  a  square,  the  pseudo  inverse  B+  is  given  by 
B+  =  (BtB)-1Bt.  Accordingly,  the  approximated  solution  of  u(t)  is 
calculated  by 


x  =  f(x,  t)  +  B(x,  t)u(t)  +  d(t ) 
y  =  [o  l  }x 

where 


0 

1 

_C 


1 

I 

1  1 

C  To 


[*i  *2]r; 


1  T 


LW  ° 


(25) 


(26) 


u(t)  =  B+(x(t),  t)[-/(x(t),  t)  -  h(x(t),t)  -d(t)+xm(t)  -Aee(t)\ 

(22) 

Then,  the  control  input  based  on  TDC  law  can  be  calculated  by 

u(t)  =  B  (x,  t)[  t)  -  k(t  -  T )  +/(x(t  -T),t  -  T ) 

+  B(x(t  -  T),  t  -  T)u(t  -  T)  +  km(t)  -  Aee(t)} 

(23) 

More  detailed  expressions  for  TDC  can  be  referred  from  Refs. 
[32,33].  To  apply  TDC  algorithm  to  PEMFC-based  DC/DC  converters, 
the  simplified  models  can  be  used  without  sacrificing  the  control 
accuracy.  Assuming  that  the  electrical  elements  are  all  ideal  and 
very  small,  i.e.,  ESR  RL  =  Rc  =  0,  the  state  variables  are  changed  as 


<*i(b 


-  w  0 


d2(t) 


[0  0]r 


The  subscript  1  represents  boost  converter  and  2  denotes  a  buck 
converter  for  the  matrices  of/,  B ,  and  d. 

Then,  the  following  reference  variables  are  chosen: 


xm\  —  *FC,ref 
xm2  —  ^o,ref 


P2  f 
o,ref 

r0'PFC 


(27) 


The  derivative  of  the  reference  state  variable  can  simply  be 
denoted  as 


x  =  [Ipc  V0}T  =  [Xi  x2]7 


(24) 


Xm(t)  =  0 


(28) 


Hence,  the  tracking  error  and  error  dynamics  are  given  by 


Fig.  4.  Schematic  circuit  of  DC/DC  converters:  (a)  boost  converter;  (b)  buck  converter. 


Table  2 

List  of  parameters  for  DC/DC  boost  converter. 
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Parameter 

Parameter  description 

Value 

C 

Capacitance 

200  pF 

L 

Inductance 

500  pH 

Rc 

ESR  of  capacitor 

30  mQ 

Rl 

ESR  of  inductor 

30  mQ 

PWM  frequency 

PWM  frequency 

30  kHz 

Table  3 

List  of  parameters  for  DC/DC  buck  converter. 

Parameter 

Parameter  description 

Value 

C 

Capacitance 

1000  pF 

L 

Inductance 

1000  pH 

Rc 

ESR  of  capacitor 

30  mQ 

Rl 

ESR  of  inductor 

30  mQ 

PWM  frequency 

PWM  frequency 

30  kHz 

e(t)  =  xm(t)  -x(t)  (29a) 

e(t)  =  K'[xm(t)  -x(t)]  (29b) 

A  canonical  form  was  presented  to  model  the  TDC  reference 
model  in  Ref.  [37].  The  output  voltage  is  the  main  control  purpose. 
Moreover,  this  study  only  used  voltage  terms  to  construct  the 
controller.  Therefore,  the  control  input  of  the  boost  converter  is 
obtained  as 


u(t)  = 


MO 


^9jFc(t_T)_M^+^nu(t_T)-Aee(t) 


Vo(t)-UFC(t-T)+Vo(t-T)[u(t-T)-'l]-LKlv0ref-vo(t)] 


MO 


The  control  input  of  the  buck  converter  is  calculated  by 


(30) 


Fig.  5.  Block  diagram  of  TDC  for  PEMFC-based  DC/DC  converters:  (a)  with  boost  converter;  (b)  with  buck  converter. 
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U(t) 


L 

yFc(0 


Mt} 

L 


-  *Vc(f  -  T) 


Voit-T) 

L 


v^l  —  u(t  -  T)  -Aee(t) 


Vo (0  -  LiFC(t  -T)  +  vpc(t  -  T)u(t  -T)-  v0(t  -  T)  -  LI<[v0 ref  -  i>0(t)] 

t'Fc(f) 


(31) 


where  is  the  numerical  derivative  of  PEMFC  current,  and  I<  is  a 
negative  constant  to  ensure  that  error  dynamics  are  convergent.  In 
boost  and  buck  converters,  I<  equals  to  -0.8  and  -0.1,  respectively. 
To  select  the  suitable  time  delay  T,  which  usually  has  same  time 
period  as  the  sampling  time  Ts  or  its  integer  multiple,  a  simple 
polynomial  characteristic  equation  is  derived  using  a  second-order 
Pade  approximation.  Time  delay  T  is  chosen  as  T  =  Ts  =  1  ms  in  our 
study.  The  TDC  law-governed  PEMFC-based  DC/DC  converter  con¬ 
trol  scheme  is  shown  in  Fig.  5. 

3.2.  Conventional  design  based  on  PIC  and  MPC 

To  test  the  superiority  of  TDC  performance,  PIC  and  MPC  are 
chosen  with  the  same  control  specification  utilized  in  the  TDC 
design.  The  transient  specifications  are  as  follows:  response  time, 
<0.2  s;  and  small  percentage  of  overshoot.  Both  PIC  and  MPC  are 
designed  based  on  the  small-signal  averaging  models  obtained 
from  the  linearization  of  Eqs.  (15)— (17).  For  the  boost  converter  [19] 


Vo(S)  =  Vp(l  -  Up)  (1  fe(1-{/q)2  5)  0  +  ^  +  RcCS) 

U(s)  LC  'rl+Rc(1-U,,)+R,,(1-Uo)2  .  RqC(Rl+Rc(1-Uq))+L5  s2 

R0LC  RotC 

(32) 

where  V0  and  U0  are  the  nominal  values  of  voltage  and  PWM  duty 
ratio,  respectively.  The  nominal  operation  point  for  boost  converter  is 
chosen  by  considering  PEMFC-rated  power  range  of  approximately 
( 14.6  V,  6  A),  which  is  given  by  VFc  =  14.6  V,  V0  =  24  V,  and  U0  =  0.392. 
The  transfer  function  at  the  nominal  operation  is  obtained  as 


V0(s)  _  -0.1801s2  -  2.9279  x  104s+  1.4659  x  10s 
0 (S)  "  s2  +  8.5694  X  102s  +  3.7700  x  106 

For  the  buck  converter  [19] 


(33) 


Frequency  (l-fe) 


Fig.  6.  Frequency  responses  of  PIC  compensated  systems:  (a)  boost  converter;  (b)  buck  converter. 


Fig.  7.  Block  diagram  of  PIC  for  PEMFC-based  DC/DC  converters:  (a)  with  boost  converter;  (b)  with  buck  converter. 
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Table  4 

List  of  MPC  parameters. 


Parameter 

Parameter  description 

Boost  DC/DC 
converter 

Buck  DC/DC  converter 

P 

Prediction  horizon 

10 

14 

M 

Control  horizon 

1 

2 

(p 

Output  weighting 

1 

0.5 

T 

Input  weighting 

1 

1 

V0,min 

Minimum  ripple  voltage 
at  nominal 

18  V 

10  V 

lFC,min 

Minimum  ripple  current 
at  nominal 

4.5  A 

2  A 

Vo, max 

Maximum  ripple  voltage 
at  nominal 

30  V 

15  V 

he, max 

Maximum  ripple  current 
at  nominal 

12.5  A 

14  A 

A^min 

Minimum  control  action 

-0.5 

-0.5 

Alimax 

Maximum  control  action 

0.5 

0.5 

V0(S)  _  Vo 


(1  +RcCs) 


°(s)  Uo  1  +  (wTRo- + RcC  +  c)s  +  |±|LCs2 


(34) 


with  the  nominal  operation  point  chosen  as  I/pc  =  14.6  V,  V0  =  12  V, 
and  H0  =  0.822.  The  transfer  function  at  the  nominal  operation  is 


V0(s)  _  4.3796  x  10-4s  +  14.5985 

U(s)  _  1.0  x  10-6s2  + 5.5217  x  10-45  +  l 


(35) 


3.2.1.  PIC  design 

In  conjunction  with  test  results,  a  suitable  PIC  is  designed  for  the 
boost  converter 

Gc(s)  =  0.0083  +  ^2,  (36) 

and  for  buck  converter 

Cc(s)  =  0.11  +-y~-  (37) 

The  bode  plots  for  the  PIC-compensated  systems  are  shown  in 
Fig.  6.  For  the  boost  converter,  the  phase  margin  is  29.5°  at  the  gain 
crossover  frequency  of  approximately  373  Hz.  By  contrast,  the 
phase  margin  is  45.8°  at  gain  crossover  frequency  of  approximately 
309  Hz  for  the  buck  converters.  The  control  diagram  is  shown  in 
Fig.  7. 


3.2.2.  MPC  design 

The  main  process  of  MPC  involves  the  calculation  of  a  series  of 
control  actions  of  PWM  duty  ratio  based  on  DC/DC  converter  model 
predictions  to  minimize  the  error  between  output  voltage  and 
setpoint  value.  The  obtained  control  sequence  calculates  new  pre¬ 
dictions  and  further  updates  the  manipulated  variables.  The  num¬ 
ber  of  new  predictions  is  called  prediction  horizon  P,  and  the 
updated  control  sequence  number  is  called  control  horizon  M, 
where  P>M.  The  linearized  model  plants  of  DC/DC  converters  with 
controlled  auto-regressive  integrated  moving  average  predictor 
(CARIMA)  [40]  are  used  to  predict  the  future  outputs.  CARIMA  is 
given  by 


a(z  1) -y(k)  =b(z  1 yu(k-\)+D(z  1 )-w(fe)-Hg-c(z  ^-e(k) 

(38) 

where  y(k),  u{k  -  1 ),  w(/<),  and  e(/<),  are  the  plant  output  prediction, 
manipulated  variable,  system  disturbance,  and  noise  at  time  instant 
k,  respectively.  The  operator  A  is  defined  by  A  =  1  -  z_1.  Moreover, 
through  model  predication,  the  manipulated  variable  is  the  solu¬ 
tion  of  the  constrained  optimization  problem  and  subjects  the  cost 
function  formulated  by 

P  n  M  n 

J  =  E  +j\k )  -  r(k  +  j)}2  +  E  Y ,-[u(k  +j  -  1)] 

ymin  <  9(k  +  j\k)  <  ymax 
^min  <  u(fe  +j\k)  <  Wmax 
Aumin  <  Au(k  +j\k)  <  A u  max 

Quadratic  programming  problem  can  be  solved  by  minimizing 
the  cost  function 

u  =  argmin J(u,r,x)  (40) 

where  y(k  +j\k)  -  r(k  +  j)  denotes  the  output  error  of  predicted 
time  instant  k  +  j  based  on  the  s  plant  information  of  time  instant  k, 
A u  is  the  control  action,  and  u*  is  the  optimal  control  input.  0j  and  jj 
represent  output  weight  factor  and  manipulated  variable  weigh 
factor,  respectively.  The  constraints  in  Equation  (39)  are  made  ac¬ 
cording  to  the  physical  limitations  and  performance  specification, 
for  example,  the  maximum  operation  voltage  of  capacitor  and 
maximum  ripple  voltage  used  in  the  boost  converter.  The  MPC 
parameters  are  obtained  through  heuristics  and  are  given  in 
Table  4.  Fig.  8  shows  the  disturbance  rejection  performance  of  MPC 
under  the  prediction  horizon  P  and  the  control  horizon  M.  The 
horizons  can  be  found  using  the  closed-loop  prediction  of  the 


Time  (9ec) 


Time  (sec) 


Fig.  8.  Prediction  horizon  P  and  control  horizon  M  tuning  under  step  load  injection  for  MPC:  (a)  boost  converter;  (b)  buck  converter. 
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Fig.  10.  Steady-state  input  condition  results  of  PEMFC-based  DC/DC  converters  without  external  applied  load:  (a)  boost  converter  case;  (b)  buck  converter  case. 


output  and  its  sensitivity  [41  ].  The  overall  MPC  control  system  of 
the  converter  is  exhibited  in  Fig.  9. 

4.  Numerical  simulation  results  and  experimental  validation 

The  numerical  simulation  results  and  experimental  validations 
are  presented  in  this  section,  based  on  the  design  parameters  in 
Tables  2  and  3.  The  main  purpose  of  the  present  study  is  to 
design  and  test  the  proposed  control  techniques  for  PEMFC-fed 
DC/DC  converters.  Tracking  24  V  for  boost  converter  and  12  V 
for  buck  converter  are  chosen  as  the  control  specifications  for  the 
simulation  and  experimental  implementation.  The  maximum 
working  power  of  100  W  for  the  converter  is  chosen  to  avoid 
concentration  loss  of  PEMFC  and  to  protect  the  PEMFC.  Simula¬ 
tions  are  performed  first  under  MATLAB/Simulink  platform, 
which  depends  on  the  designed  control  schematic  diagrams 
shown  in  Figs.  5,  7  and  9.  Real-time  experiments  were  also  per¬ 
formed  with  NI  DAQ  system  with  LabVIEW  to  validate  the 
effectiveness  of  the  proposed  control. 

4.1.  Results  of  numerical  simulation 

The  results  of  the  three  control  methods  under  steady-state 
input  conditions  without  externally  applied  load  are  shown  in 
Fig.  10.  Each  controller  shows  desired  control  performance,  i.e.,  the 
settling  time  is  less  than  0.2  s,  and  tracking  error  is  within  the 
specified  level.  Moreover,  TDC  shows  the  best  control  performance 
with  the  fastest  response  and  without  overshoot  in  both  boost  and 
buck  converters,  whereas  PIC  shows  sluggish  response,  and  MPC 
exhibits  large  overshoot.  An  abrupt  external  load  change  is  also 
applied  to  determine  the  disturbance  rejection  ability.  Even  in  the 
presence  of  an  external  disturbance,  DC/DC  converters  are  still 
required  to  supply  the  constant  voltage  accurately  and  promptly  to 
the  secondary  power  buses.  The  results  of  the  boost  converter 
under  abrupt  external  disturbance  engagement  are  shown  in  Fig.  11. 
TDC  shows  the  most  robust  tracking  performance  against  the  other 
controllers.  Immediately  after  the  introduction  of  the  disturbance, 
TDC  vanishes  DC/DC  output  voltage  ripple.  Hence,  error  diminishes 
rapidly.  For  example,  when  the  maximum  disturbance  occurs,  TDC 
settling  time  is  just  approximately  0.01  s  with  4.6  V  overshoot. 
However,  MPC  settling  time  is  approximately  0.04  s  with  10.8  V 
overshoot  and  that  of  PIC  is  approximately  0.2  s  with  9.6  V  over¬ 
shoot.  Simultaneously,  PEMFC  shows  the  most  stable  stack  voltage 
response  under  TDC,  which  guarantees  a  stable  fuel  cell  operation. 

TDC  shows  the  best  control  performance  (Fig.  12(c)),  among  the 
three  controllers  for  PEMFC-fed  DC/DC  buck  converter  in  terms  of 
response  time,  overshoot,  and  ripples.  MPC  shows  faster  response 


in  steady-state  operation.  However,  MPC  has  higher  overshoot  and 
larger  ripples  in  the  presence  of  a  disturbance.  The  ripples  may 
impact  the  PEMFC  stack  voltage  oscillations  as  shown  in  Fig.  12(b), 
which  may  degrade  the  fuel  cell  performance  in  long-term  usage. 
By  contrast,  the  disadvantage  of  PIC  is  its  sluggish  response  in  the 
presence  of  a  disturbance.  For  example,  at  maximum  applied  load, 
PIC  requires  approximately  0.2  s  to  track  the  setpoint.  The  simu¬ 
lation  results  reveal  that  TDC  shows  the  most  excellent  transient 
responses  for  PEMFC-fed  boost  and  buck  DC/DC  converter 
application. 


4.2.  Experimental  validation 

An  experimental  test  rig  is  constructed  using  NI  PCI-6229  DAQ 
board,  SCXI-1313A  signal  conditioner,  H-100  PEMFC,  DC/DC  boost 
and  buck  converters,  electronic  load,  and  TDS3034B  series  digital 
phosphor  oscilloscope  (Fig.  13).  PCI-6229,  a  data  acquisition  device 
with  32  analog  inputs  (16  bits,  250  kS  s-1)  and  4  analog  outputs 
(16  bits,  833  kS  s-1),  provides  optimized  functionality  for  cost- 
sensitive  applications  [42].  SCXI-1313A  is  used  for  high-voltage 
attenuation  and  AC/DC  coupling.  NI  LabVIEW  system  provides  a 
widely  used  graphic  user  design  interface,  which  is  easy  to  program 
with  various  and  friendly  hardware  I/O  ports.  Therefore,  real-time 
measurement  and  control  for  PEMFC  converter  systems  are  easily 
implemented.  The  voltage  and  current  sensing  are  recorded 
through  oscilloscope  and  LabVIEW  system.  The  current  is  trans¬ 
formed  to  voltage  through  a  Hall-effect  sensor,  WCS2705,  or 
WCS2720. 

Fig.  14  shows  the  response  of  the  starting  period  of  the  PEMFC- 
fed  DC/DC  boost  converter,  where  the  five  fluctuated  PEMFC 
voltage  loops  are  caused  by  H-100  PEMFC  built-in  start  mode.  H- 
100  PEMFC  has  no  relation  with  the  DC/DC  converter  control  effect. 
Large  fluctuations  caused  by  PEMFC  “short  circuit”  operation  mode 
at  the  start  period  that  may  increase  stack  voltage  [38].  At  this  time, 
TDC  shows  a  rapid  response  (<0.1  s)  with  smaller  overshoots 
(maximum  overshoot  of  3  V)  to  track  the  setpoint,  whereas  MPC 
shows  the  worst  performance  with  large  overshoots  of  approxi¬ 
mately  7  V. 

The  disturbance  rejection  ability  was  tested  through  abrupt 
external  load  engagement.  The  results  of  the  boost  converter 
controlled  by  PIC,  MPC,  and  TDC  are  shown  in  Fig.  15(a)— (c), 
respectively.  TDC  exhibits  the  best  transient  performance  in 
considering  the  response  time  and  overshoot,  whereas  MPC 
shows  corresponding  values  of  0.2  s  and  4.1  V,  and  those  of  PIC 
are  0.3  s  and  7  V  under  maximal  load.  Fig.  16(a)— (c)  shows  the 
response  of  disturbance  rejection  ability  for  the  DC/DC  buck 
converter  under  sudden  external  load  change.  TDC  displays  the 
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Fig.  11.  Response  of  step  load  applied  to  PEMFC-based  DC/DC  boost  converter:  (a)  with 
PIC;  (b)  with  MPC;  (c)  with  TDC. 


best  performance  against  PIC  and  MPC.  Considering  the  results 
presented  in  Figs.  15  and  16,  TDC  is  the  most  suitable  control 
technique  among  the  three  controllers  for  PEMFC-fed  DC/DC 
converters.  Moreover,  a  ramped  load  is  applied  to  the  converters 
to  examine  the  operation  range  of  the  converter.  The  tests  are 
implemented  on  the  boost  converter  using  PIC  and  TDC  only.  PIC 
and  TDC  almost  showed  similar  control  results  (Fig.  17).  Both 
have  quick  response  times  and  small  ripples,  which  indicate  that 
both  controllers  are  useful  in  regulating  the  output  voltage  for 
the  variation  in  ramp  load.  TDC  result  shows  much  faster 


transient  dynamics  for  boost  converter;  however,  better  transient 
response  cannot  be  obtained  for  buck  converter  even  TDC  is 
applied.  This  is  because  boost  converter  has  the  non-minimum 
phase  characteristics. 

The  tracking  errors  for  start  mode  and  step-load  change  oper¬ 
ation  mode  are  calculated  by  root  mean  square  error  (RMSE).  RMSE 
comparisons  among  PIC,  MPC,  and  TDC  are  listed  in  Table  5.  TDC 
shows  the  minimum  errors  to  track  the  setpoint  in  boost  and  buck 
converters  under  different  working  conditions. 


Fig.  12.  Response  of  step  load  applied  to  PEMFC-based  DC/DC  buck  converter:  (a)  with 
PIC;  (b)  with  MPC;  (c)  with  TDC. 


Y.-X.  Wang  et  al.  /  Journal  of  Power  Sources  261  (2014)  292-305 


303 


Fig.  13.  Picture  of  an  experimental  test  rig. 


Fig.  14.  Experimental  results  of  PEMFC-based  DC/DC  boost  converter  with  start  operation:  (a)  with  PIC;  (b)  with  MPC;  (c)  with  TDC.  Notes  on  the  oscilloscope:  Channel  1  is  the 
PEMFC  stack  voltage;  Channel  2  is  the  DC/DC  converter  output  voltage;  Channel  3  is  the  transformed  voltage  of  PEMFC  stack  current;  Channel  4  is  the  transformed  voltage  of  load 
current. 


Fig.  15.  Experimental  results  of  PEMFC-based  DC/DC  boost  converter  step  load  injection:  (a)  with  PIC;  (b)  with  MPC;  (c)  with  TDC.  Notes  on  the  oscilloscope:  Channel  1  is  the  PEMFC 
stack  voltage;  Channel  2  is  the  DC/DC  converter  output  voltage;  Channel  3  is  the  transformed  voltage  of  PEMFC  stack  current;  Channel  4  is  the  transformed  voltage  of  load  current. 


5.  Conclusions 

This  study  introduces  a  robust  control  technique  for  PWM- 
regulated  DC/DC  converters  used  for  PEMFC  application.  The 
study  focuses  on  the  PEMFC-fed  PWM  DC/DC  converter  controller 
to  supply  both  constantly  high  (boost  action)  and  low  (buck  action) 
DC  voltage  to  the  subsequent  power  link.  PEMFC  dynamics  is 
modeled  through  equivalent  electrical  circuits  with  a  variable 


capacitor,  and  linearized  DC/DC  converter  models  are  added  for  full 
control.  Based  on  the  control-oriented  full  model,  including  PEMFC 
and  two  DC/DC  converters,  three  control  techniques,  namely,  TDC, 
MPC,  and  PIC,  are  designed  and  simulated  using  MATLAB/Simulink. 
The  simulation  results  reveal  that  TDC  is  the  most  preferable  con¬ 
trol  technique  for  boost  and  buck  DC/DC  converters  in  PEMFC 
application.  TDC  provided  a  stable  and  fast  DC  link  voltage.  Boost 
and  buck  DC/DC  converters  are  constructed  and  connected  to  the 
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Fig.  16.  Experimental  results  of  PEMFC-based  DC/DC  buck  converter  step  load  injection:  (a)  with  PIC;  (b)  with  MPC;  (c)  with  TDC.  Notes  on  the  oscilloscope:  Channel  1  is  the  PEMFC 
stack  voltage:  Channel  2  is  the  DC/DC  converter  output  voltage;  Channel  3  is  the  transformed  voltage  of  PEMFC  stack  current;  Channel  4  is  the  transformed  voltage  of  load  current. 


Fig.  17.  Experimental  results  of  PEMFC-based  DC/DC  boost  converter  ramp  load  injection:  (a)  with  PIC;  (b)  with  TDC.  Notes  on  the  oscilloscope:  Channel  1  is  the  PEMFC  stack 
voltage;  Channel  2  is  the  DC/DC  converter  output  voltage;  Channel  3  is  the  transformed  voltage  of  PEMFC  stack  current;  Channel  4  is  the  transformed  voltage  of  load  current. 


Table  5 

Results  of  RMSE  with  experimental  results. 


Operation  mode 

PIC 

MPC 

TDC 

Start  mode  with  boost  converter 

7.5825 

7.9273 

7.5977 

Start  mode  with  buck  converter 

5.4204 

5.3570 

5.0889 

Step  load  mode  with  boost  converter 

0.9733 

0.4357 

0.2844 

Step  load  mode  with  buck  converter 

0.3195 

0.4840 

0.1384 

H-100  PEMFC  to  prove  the  efficacy  of  the  TDC  experimentally.  NI 
DAQ  LabVIEW  system  is  utilized  as  a  real-time  control,  and  a  series 
of  the  tests  are  implemented  to  compare  the  effects  of  the  three 
controls.  The  experimental  results  showed  the  superiority  of  TDC 
performance  to  those  of  PIC  and  MPC  in  providing  stable  and  fast 
DC  power  link,  even  when  the  input  voltage  or  external  output  load 
varies.  Its  relatively  easy  design  process  is  another  advantage  of 
using  TDC  because  tuning  its  parameters  is  much  simpler  than 
those  of  PIC  and  MPC.  The  proposed  topology  and  control  approach 
can  be  utilized  for  other  power  source-based  fed  converters,  such  as 
photovoltaic  cells,  super  capacitors,  and  batteries. 
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